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Q ABSTRACT 



The prompt emission of gamma-ray bursts (GRBs) is widely thought to be 
radiation from accelerated electrons, but an appreciably larger amount of energy 
could be carried by accelerated protons, particularly if GRBs are the sources of 
ultra-high-energy cosmic rays (UHECRs). We model the expected photon spectra 
+-> ■ for such "proton-dominated" GRBs in the internal shock scenario through Monte 

Carlo simulations, accounting for various processes related to high-energy elec- 
trons and protons. Besides proton and muon synchrotron components, emission 
£> . from photomeson-induced secondary pair cascades becomes crucial, generally en- 

hancing the GeV-TeV and/or eV-keV photons and offering a signature of UHE 
Q\ . protons. In some cases, it can overwhelm the primary electron component and 

^ . result in GRBs peaking in the 10 MeV - 1 GeV range, which may be relevant 

to some bursts discussed in a recent re-analysis of EGRET TASC data. The 
dependence of the spectra on key quantities such as the bulk Lorentz factor, 
magnetic field and proton-to-electron ratio is nontrivial due to the nonlinear na- 



% 



ture of cascading and the interplay of electron- and proton-induced components. 
Observations by Fermi, ground-based telescopes and other facilities should test 
these expectations and provide critical constraints on the proton acceleration 
efficiency. 

Subject headings: cosmic rays — gamma rays: bursts — gamma rays: theory - 
radiation mechanisms: nonthermal 
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Introduction 



The prompt emission of gamma-ray bursts (GRBs) is believed to arise f rom ultrarel- 



ativistic outflow s with bulk Lorentz factors T > 100 (see, e.g., reviews by iPiran 1 12005 



Meszaros I I2006l l 



In the popular internal shock model, collisions among inhomogeneities 
within the flow lead to formation of shocks that convert a fraction of the bulk kinetic energy 
into Fermi-accelerated r elativistic electrons, who se synchrotron emission powers the observed 
MeV-band gamma-rays (IRees fc Meszaros II 19941 ). Initially, most of the kinetic energy as well 
as the internal energy generated via shock dissipation are likely carried by protons, so such 
models entail the operation of a physical mechanism that transfers energy from protons 
to electrons on sufficiently short timescales. This presumably occurs via collective electro- 
magnetic processes, as simple Coulomb collisions may be too slow. A general problem in 
collisionless shock theory and GRB models in particular is that this mechanism is poorly 
understood, and one must frequently resort to a phenomenological parametrization. In view 
of the large observed energy in MeV gamma rays, the efficiency of proton-to-electron energy 
transfer is usually considered to be high. However, this is by no means physically guaranteed. 
In the case of supernova remnant shocks, the total energy in acce lerated electrons is ofte n 
constrained observationally to be much less than in protons (e.g. lAharonian et al. 1120061 ). 
Since we do not yet understand the nature and total energy budget of the central engine, 
we cannot readily exclude the possibility that GRBs actually contain a significantly larger 
amount of energy in protons compared to that radiated by the accelerated electrons. 

Furthermore, a natural expectation is that the shocked protons are also Fermi-accelerated 
The physical conditions in internal shocks may allow maximum energies > 10 20 eV , so GRBs 



are potential sources of the observed ultra -high-energy cosmic rays (UHECRs; IWaxman 



19951 ; IVietri I Il995l ; iMilgrom fc Usov 1119961 ). The total energy in accelerated protons that 
must be suppl i ed pe r burst depends on a number of uncertain factors (see also App. B of 
Murase et al. I (120081 )). The required local UHECR emissivity at prot o n energy e p ~ 10 19 eV 
is e 2 p dN p /de p ~ 0.8 x 10 44 erg Mpc^yr" 1 dWaxman k Bahcall I fl998J : bermer 1 120071 1 . Post- 



S WIFT estimates of the local rate of long GRBs range from 0.2 — 1 Gpc 3 yr l if the GRB 
rate is proportional to the star formation rate, down to ~ 0.05 Gpc~ 3 yr~ 1 if the GRB rate 



evolves more strongly with redshift, which may be observationally favored (e.g. iDaigne et al 



20061 ; iLe fc Dermer 1 120071 : iGuetta Sz Piran 1120071 ). Assuming a power-law proton spectrum 
with index p p = 2, the necessary isotropic-equivalent energy per burst in accelerated protons 
integrated over e p ~ 10 9 — 10 20 eV is E p ~ 2 x 10 54 — 3 x 10 55 erg, which is approximately inde- 
pendent of the actual beaming factor. Steeper spectra and h ence even larger E„ are called for 
if GRBs also contribute significantly to CRs below 10 19 eV ( jWick et al. II2004I ). To be com- 
pared is the corresponding energy in accelerated electrons E e , which can be roughly equated 
with the observed, isotropic-equivalent MeV gamma-ray energy Kyi S0 , typically ~ 10 53 erg 
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and up to ~ 10 54 erg in the 1 — 10 4 keV rest-frame band fjKocevski Sz Butler 1120081 ). Thus, 
in order for GRBs to be viable sources of UHECRs, the latest observations point to a highly 
proton-dominated energy budget, E p /E c > 10 — 100. The observed heterogeneity of GRBs 
also suggests that not all bursts may be equally efficient UHECR accelerators, in which case 
even higher E p /E e may be warranted for a subset of the bursts. 

It is therefore of great interest whether such "proton-dominated" GRBs can be di- 
agnosed observationally. A promising window is GeV-TeV gamma-rays, where distinctive 
signatures of UHE proton acceleration may show up, such as synchrotron emiss ion from pro 



tons, muons or secondary particles injected via photomeson interaction s (e.g. IVietri 1119971 ; 



Bottcher fc Dermer 1119981 : iGupta fc Zhang 1 120071 : lAsano fc Inoue 1 120071 . hereafter AI07). 



AI07 recently undertook a detailed investigation of such emission processes utilizing a 
comprehensive Monte Carlo code. However, having assumed that the accelerated protons 
do not carry excessive extra energy, their study was restricted to E p /E c = 1. In view of 
the above possibilities, here we follow and extend the work of AI07 to E p /E c > 1. The 
results, which are often qualitatively and drastically different from AI07, are discussed in 
relation to existing and upcoming observations. Note that high-energy emissio n from proton 
dominated GRBs has b een discussed previously in different contexts (e.g. 
Asano fc Takaharal 120031 ) . 



Totani 1998 



After a recap of our formulation in §2J we discuss the results and their observational 
implications in §3] and §H respectively, and conclude in §5J 



Model and Methods 



We briefly summarize the model and methods of AI07, which should be consulted for 
more details. In accord with the internal shock paradigm, the emitting region corresponding 
to an individual pulse in the prompt light curve is taken to be a homogeneous shell expanding 
with T at radii R from the central engine. The comoving width of the shell is I = R/T and 
the pulse timescal e in the observer fra me is At = R/T 2 c, as long as R exceeds the shell 
spreading radius (IMeszaros et al.lll993l ). which is always the case here. Shock dynamics 
and time variability are not explicitly treated, so our results should be interpreted as the 
time- averaged spectra for each pulse. 

With given injection of accelerated electrons and protons in magnetic field B, we solve 
self-consistently for the distribution of particles and photons in the shell using Monte Carlo 
techniques. The ti me steps are always taken to be sufficiently shorter than the particle 
cooling timescales (lAsano 1 120051 ). In addition to synchrotron and inverse Compton (IC) 
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emission from all particles, our code includes synchrotron self-absorption, cascade processes 
with photon-photon (77) production of electron-positron pairs (e 1 * 1 ) and Klein-Nishina regime 
Compton scattering, as well as proton-induced processes such as photomeson (jry) interac- 
tions and secondary pion, muon and pair injection. We adopt experimental results for the 
cross sections of pry — > mr + , ptr , nir + 7r° and p7r + 7r~, while p~y — > pir 7r° is neglected in view of 
its small cross section. In case the primary proton is converted to a neutron, we assume that 
it continues to interact with photons in the shell during the comoving expansion timescale 
texp = l/c. We do not accco unt for the minor contribution from neutron-decay electrons 
( jRazzaque fc Meszaros 1120061 ). More de t ails on the treatment of meson p roduc tion and their 
decay products can be found in lAsano I (120051 ) and lAsano fc Nagataki I (120061 ). 



Furthermore, we now account for the Bethe-Heitler (BH ) pair production pro c ess (jr y 



pe + e ), whose cross section and inelasticity are taken from IChodorowski et al. I (Il992l ). In 



the present context, the proton energy loss is always dominated by photopion production, and 
the huge compactness of GRBs implies that the resultant electromagnetic cascade emission 
is not very sensitive to the details of particle injection at high energies. Thus, compared 
to cases neglecting the BH process, we find that its inclusion here only leads to modest 
enhancements of the secondary photon emission, by at most a few tens of percent. 

Primary electrons with total energy density U e are injected with a power-law distribution 
n e (7e) oc 7~ Po in the range of Lorentz factors 7 c , m in < 7 C < 7 ,max- The balance of Fermi 
acceleration and radiative cooling timescales gives 7 e , m ax- Likewise, protons with total energy 
density U p are injected with a distribution n p (7 p ) ex 7 P Pp in the range 7 Pim i n < 7 P < 7 Pi max- 
We obtain 7 Pimax by equating t acc = 'j p m p c 2 /eBc, the Fermi acceleration time in relativistic 
shocks, to minftexp, tj oss), where t\ oss is the energy loss time due to synchrotron, IC and pj 
cooling (lAsano 1120051 ). In mildly relativistic internal shocks, 7 Pjm i n should be of order unity; 
here we take 7 Pim i n = 10. 

The injection index for electrons is fiducially chosen to be p e = 2.5, implying (3 ~ 2.25 
for the spectral index above the synchrotron peak energy. This is consistent with the mean 
of the (3 values measured by BATSE, albeit with a considerable dispersion, from (3 < 1.5 to 



(3 > 3.0 (IPreece et al. Il2000l ; iKaneko et al. 1120061 ). For protons, our fidu cial index is p p = 2.0, 



appro priate when GRBs contribute to UHECRs only above 10 eV ( jWaxman fc Bahcall 



19981 ); steeper spectra would increase still the energy demands. Note that the values of p e 
and p p relevant to our results each correspond to very different energy ranges; GeV-TeV 
for electrons and 10-100 PeV for protons in the comoving frame. Although the injection 
spectra for the two species are expected to be the same at low energies where their gyroradii 
overlap, p c > p p may be effectively realized if the proton spectrum covering 7-8 decades in 
energy deviates from a pure power-law and becomes concave. This may plausibly occur due 
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to 1) nontrivial geom etry and wavelength distribution of magnetic turbulence at the shock 



(Niemiec et al. 



( jBaring fc Kirk 



1991 



2006 ) , 2) nonlinear back-rea ction of CR pressure on the shock structure 



Malkov fc Drury Il200ll ). or 3) superpositio n of pre-existing and newly - 



injected particles originating from different regions in the outflow fjBosnjak fc Daigne 1120081 ). 
Nevertheless, in view of the observed spread in j3 and the uncertainties associated with 
obtaining time-integrated spectra, we also discuss cases with p c = p p in §3.21 



Some combinations of the remaining parameters are constrained so as to reproduce 
typically observed properties of the MeV emission. For given B and T, 7 eim i n is chosen such 
that the observed synchrotron peak energy for nearby bursts is e p k = Tj^ min heB/m e c ~ 
300 keV. Since the fast-cooling, primary electrons radiate away most of their energy as 
MeV photons within At, E e = (4irT 2 R 2 cAt)U c ~ {4irR 3 )U c can be identified with E sh , the 
observable, isotropic-equivalent MeV pulse energy. 

Instead of U c , U p and Ub = B 2 /8n, hereafter we use e e , e p and ee, the conventional 
parametriza tion of the corre sponding energies as fractions of the shock-dissipated internal 



energy (e.g. 



Meszaros 1 120061 ) . Thus e B /e = Us/U e and e p /e c = U p /U c = E p /E e . In place 
of R, we choose the observable At as a parameter and set At = 0.1 s for simplicity. Below 
we only show the spectra corresponding to single pulses. For bursts composed of N similar 
pulses, the duration-integrated energy would be simply N times larger, E % i so = NE^. The 
set of parameters are then At, E s ^, T, ee/ce and e p /e e . All spectra are plotted as observed 
fluence versus photon energy, assuming a GRB redshift z = 0.1. Spectral attenuation by 
intergalactic 77 absorption is neglected. 



Results 



3.1. Fiducial Spectral Indices 

First we discuss different cases with our fiducial values of p e = 2.5 and p p = 2.0. As 
mentioned above, e p k is chosen to have the typically observed value of 300 keV. Prompt 
emission spectra of single pulses for E sh = 10 51 erg, T = 300, e B /e c = 1 an d varying 
e p/ e c = 10 — 100 are shown in Figured) The sharp spectral cutoffs at low and high energies 
are due to synchrotron self-absorption and 77 absorption, respectively. This applies to 
all spectra below when such sharp cutoffs are seen. Most remarkable is the prominent e^ 
cascade component, i.e. synchrotron and IC emission from secondary e ± triggered by p-f 
interactions of UHE protons with low energy photons. For the lower range of e p /e e , primary 
synchrotron photons constitute the main jry target. However, as the proton content increases, 
the target photons become dominated by synchrotron emission from the low energy part of 
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the secondary e^ themselves. The dependence of the spectra on e p /e is therefore nonlinear 
and not simply proportional, as apparent in Figure [TJ The secondary photons also affect the 
primary synchrotron component (dashed curves in Figure [T]) through enhanced IC cooling, 
even though the injection distribution is unchanged. 
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Fig. 1. — Single pulse, prompt photon spectra for varying e p /e e as labeled. Other parameters 
are marked above the figure. Dashed curves denote the primary contribution only, whose 
peak flux decreases with e p . Dot-dashed curves denote separately the electron synchrotron 
(labeled eSY) and inverse Compton (elC) components without 77-absorption effects for 
e p /e e = 100. 

In general, cascade emission significantly hardens the high-energy spectra. Since sec- 
ondary e with Lorentz factors < 7 Cjm i n can be injected in the cascade, it can also give rise 
to excess UV-to-X-ray emission lying above the extrapolation of the sub-MeV spectra, as 
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seen for e p /e c = 10 — 30 in Figured! The entire spectra thus tends to become flat in ef(e). 

The case of e p /e c = 100 is drastically different. Here the proton-induced secondary 
emission totally overwhelms any primary electron component, resulting in a hard spectrum 
peaking at 10-100 MeV. Although approximately a single power-law between 100 eV and 30 
MeV, in fact it comprises two emission processes by secondary e ± , mainly synchrotron < 
MeV and IC > MeV (dot-dashed curves in Figured]). Despite e B /e e = 1, IC can dominate 
over synchrotron since the energy density of secondary e exceeds both U-q and U e . 

The comoving photon density n 7 is decisive for both 1) the 77 optical depth r 77 and 
hence the 77 cutoff energy e 77 , and 2) the efficiency of py interactions and hence the sec- 
ondary cascade emission. Figure [2] displays single pulse spectra for T = 300, e B /e c = !> 
e p /e c = 10, and varying pulse energies E s ^ = 10 49 — 10 51 erg. Higher E s ^ implies higher ri 7 , 
and consequently stronger jry components as well as lower e 77 . Since n 7 oc T~ 5 with other 
parameters fixed, varying V has larger effects. Shown in Figure [3] are single pulse spectra 
for E sh = 10 50 erg, e B /e c = 1, e p /e e = 30 and T = 100 - 1000. T = 100 allows a high e pk , 
cascade-dominated spectrum, even though e p /e e is 3 times less than the analogous case in 
Figured! Increasing T leads to higher maximum energies and less cascade contribution. The 
spectral hardening > 0.1 GeV for T = 300 and > 10 GeV for T = 1000 is due to secondary 
IC. 

Thus high proton-dominance does not always result in conspicuous proton-induced emis- 
sion if T is sufficiently high. Conversely, the absence of hard, high-energy components does 
not necessary rule out proton-dominated GRBs. In fact, the conditions most favorable for 
contributing to UHECRs is that they escape the source with minimal p~y losses, which cor- 
responds roughly to the criterion T > 300(At/0.1s) _0,3 (i? s h/10 51 erg) 0,2 in our model (AI07). 
On t he other hand, T can b e observationally constrained through its strong influence on e 77 
(e.g- " 



Lithwick fc Sari I l200ll . AI07). Since the pulse energy E s ^ and timescale At are also 
measurable, we may hope to identify bursts where pj losses are likely to be efficient, and 
then constrain e p /e e from the high-energy spectra, although some degeneracy with ee/ce will 



remain. 



Figure |H shows single pulse spectra for E sh = 10 51 erg, T = 300, e p /e c = 30 and varying 
^bAg = 0.1 — 10. The ee/ec = 1 case is the same as in Figured! Higher B causes steeper 
spectra with stronger secondary synchrotron relative to secondary IC, while lower B is vice- 
versa and produces a 100 MeV peak spectrum. However, the dependence on B can also be 
nontrivial. In Figure \5[ we show spectra for E s ^ = 10 51 erg, T = 1000, e p /e G = 100, and 
varying ee/ce = 0.1 — 100 (note that cbAp < !)• The higher T allows spectra extending 
into the TeV regime, but renders P7 processes inefficient despite the high proton-dominance. 
All cases exhibit spectral bumps around 0.1-1 TeV, but their origins are quite different. For 
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Fig. 2. — Single pulse, prompt photon spectra for varying E 8 h as labeled. Other parameters 
are marked above the figure. 



ee/ec < 1, this is due to secondary e ± IC, which is weaker for higher B. However, when 
e B,/ e c ^ 10, the bump is stronger again, owing to the appearance of synchrotron emission 
from protons and muons, their ratio being roughly 2 to 1 for e-QJe c = 10 (dot-dashed curves 
in Figure [5]). For €b/€ c = 100, we obtain a pronounced proton synchrotron TeV peak, as 
well as enhanced emission at lower energies from synchrotron radiation by e^ produced via 
77 absorption. 
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Fig. 3.— Single pulse, prompt photon spectra for varying T as labeled. Other parameters 
are marked above the figure. 



3.2. Equal Proton and Electron Indices 

We now consider situations with p e = p p , as would occur if the proton spectrum was a 
single power-law over its entire energy range. Similar to the above, Figure [6] testifies that the 
spectrum for e p /e c =30 and p c = p p =2.0 can result in a hard GRB with photon index ~ 2 up 
to 10 GeV. It is interesting to note that in such cases, the spectral shape around the MeV 
peak alone may not always reveal the correct value of p e . However, for p c = p p = 2.2, the 
fraction of UHE protons and the associated cascade emission is greatly diminished, except for 
a slight distortion of the spectrum above 100 MeV. The proton contribution becomes totally 
negligible for p c = p p = 2.5, for which neither UHECRs nor neutrinos are significantly 
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Fig. 4. — Single pulse, prompt photon spectra for varying ee/ee as labeled. Other parameters 
are marked above the figure. Dashed curves denote the primary components only, whose 
peak flux decreases with ee- 

generated at any rate. 



4. Observational Implications 

A unique property of proton-dominated GRBs is that their photon spectra can some- 
times manifest very high peak energies in the 10 MeV-1 GeV range due to pq cascade 
emission ( F igures [T1I31I4161) . This seems at variance with commonly observed values of e p k ~ 
0.1-1 MeV (IKaneko et al. 1120061 ). However, through a recent re-analysis of EGRET TASC 
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Fig. 5.— Single pulse, prompt photon spectra varying ee/ce as labeled. Other parameters 
are marked above the figure. Dot-dashed curves denote separately the electron synchrotron 
(eSY), proton synchrotron (pSY) and muon synchrotron (yLtSY) components without 77- 
absorption effects for e-Q,/e c = 10 and 100. 



data, iKaneko et al. I (120081 ) reported a GRB with apparently very high e pk > 170 MeV 



as well as a few others with significant high-energy excess (see also I Gonzalez et al. 1 120031 ) . 



Some studie s have also indicated pot ential observational biases against BATSE detections 
of high e p k fjLloyd fc Petrosian Ill999l ). At this moment, it is unclear how often such high 
£ p k bursts occur, and whether they are relevant to the proton-dominated cases discussed 
here, or simply reflect a primary synchrotron peak energy that is much higher than average 
(rather than the values we have assumed here). In any case, the existence an d nature of such 
bursts will be definitively probed through ongoing observations by Fermi (jOmodei 1120061 ) 
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Fig. 6.— Single pulse, prompt photon spectra for varying values of p c = p p as labeled. 
Other parameters are marked above the figure. Thick and thin dashed curves denote the 
primary components only, for p e = p p =2.0 and 2.2, respectively. 



and AGILE (ILongo et al. 1 120071 ). Note that it is also conceivable that some GRBs possess 
conservative proton energies, say E p ~ 10 53 erg, but with E p /E e > 1 so that the MeV 
emission is relatively weak. Even if unimportant for UHECRs (HI]), new generation satellites 
should also probe such MeV-weak bursts. 

The j?7 cascade can also induce excess low-energy emission (Figures [TlliZPl) . which do not 
seem typical of known GRBs. However, the y may be relevant for some BATSE bursts with 
soft excess components (IPreece et al. Ill996l ). or possibly a fraction of the X-ray rich GRBs 
(jSakamoto et al. 1120051 ). Fermi and AGILE observations of the accompanying high-energy 
excess will provide a test. 



TeV detections of GRBs have yet to be achieved (e.g. lAtkins et al. 1120051 ; I Albert et al 
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20071 ; iHoran et al. 1 120071 ; lAharonian et al. 1120091 ). but some of the components discussed 



above may be eventually observed by current ground-based facilities such as MAGIC (II), 
HESS (II), VERITAS, CANGAROO III, or the future projects CTA, AGIS, HAWC, etc. 
For example, MAGIC may detect the luminous proton synchrotron emission for ee/e c = 100 
in Figure [5] at 0.1 TeV beyo nd z ~ 1, assuming E 1 - m) = 10 53 erg and the latest estimates of 



intergalactic 77 absorption (jAlbert et al. 1120081 ). 



Distinguishing between primary electron IC and proton-induced emission components 
may not be easy from the spectral shape alone. However, since the synchrotron and/or 
photomeson cooling timescales for UHE protons are considerably longer than the cooling 
timescales for GeV-TeV emitting primary electrons, we can expect important differences in 
their variability properties, which should provide further observational clues. Although this 
work was limited to time-averaged pulse spectra, a desirable next step is to perform explicitly 
time-dependent calculations. 



Conclusions and Outlook 



Proton-dominated GRBs are motivated by physical considerations of particle accelera- 
tion in collisionless shocks, as well as their potential to be the origin of UHECRs. In GRB 
UHECR scenarios, the spectral index for protons at UHE must generally be harder than 
the typical indices for electrons emitting in the multi-MeV range, which may be possible de- 
pending on the physics of particle acceleration, cooling and/or shock formation, as discussed 
in §2j Characteristic emission signatures can then result, such as high peak energy bursts 
and/or excess low-energy emission from photomeson-triggered pair cascades, or luminous 
spectral bumps from proton synchrotron emission. If the indices for electrons and protons at 
the respective energies are equal, proton-related components may still be visible as long as 
the index < 2.2, but not for steeper spectra. Through detailed observations of spectra and 
variability, we may hope to disentangle the proton-induced components from the competing 
emission process of inverse Compton from primary electrons. 

Other observable consequences of p roton-dominated GRBs may include contributions 
to G alactic CRs (e .g. IWick et al. 1 12004 ) and the diffuse high-energy neutrino background 



[e.g. 



Murase 1120071 ). 



We note that if some GRBs actually emit stronger GeV-TeV components than previously 
expected as discussed here, they could play an increased role in probing high- z intergalactic 



radia tion fields (Inoue et al., in prep.) as well as intergalactic magnetic fields (jlchiki et al. 
20081 . and references therein). 
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